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Abstract 

In fabrication of nano- and quantum devices, it is sometimes critical to position individual dopants at certain sites 
precisely to obtain the specific or enhanced functionalities. With first-principles simulations, we propose a method 
for substitutional doping of individual atom at a certain position on a stepped metal surface by single-atom 
manipulation. A selected atom at the step of Al (111) surface could be extracted vertically with an Al trimer-apex 
tip, and then the dopant atom will be positioned to this site. The details of the entire process including potential 
energy curves are given, which suggests the reliability of the proposed single-atom doping method. 
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Background 

Single-atom manipulation, which was first introduced by 
Eigler et al. and realized experimentally on Ni (111) sur- 
face with a scanning tunneling microscope (STM) tip, 
provides a way to fabricate nanostructures with atomic 
precision [1-7]. Besides the STM tip, for nonconductive 
surface, the tip of an atomic force microscope (AFM) 
has also been applied to achieve various single-atom 
manipulations [8-10]. Studies show that merely by the 
mechanical interaction force acting between the tip and 
atom, complex manipulations can still be accomplished 
besides the primary lateral and vertical manipulations. 
For instance, on Al (111) surface, a reversible modifica- 
tion of the configuration of supported nanoclusters with 
atomic precision by tip was demonstrated in our previ- 
ous simulations [11]. Also, the work on Si (111) surface 
given by Sugimoto et al. shows that an atom from the 
AFM tip can interchange with a surface adatom in a re- 
versible exchange procedure [9]. Through this vertical 
manipulation, a single Si atom can be precisely posi- 
tioned into or extracted from the Sn layer. As the size of 
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devices shrinks to nanoscale or even to atomic scale, be- 
sides configuration of nanostructure, the number of iso- 
lated atoms of certain species and their location could 
modify their functionality and performance [12,13]. There- 
fore, it is sometimes demanded to position dopants at cer- 
tain sites precisely. For example, by STM-based hydrogen 
lithography, a single-atom transistor in which an individual 
P dopant atom has been placed within a silicon device 
with a spatial accuracy of one lattice site was demonstrated 
recently [14]. In another work, by Cs atom doping with a 
STM tip, spin of individual magnetic molecules as basis of 
quantum computer was successfully controlled [15]. 

On metal surfaces, influences of tip structure on the 
manipulation were intensively investigated in our previous 
work [16], and it was shown that the trimer-apex tip, a 
model of blunt tip in the experiment, is capable of trans- 
forming the configuration of the Al nanocluster reversibly 
[11]. The specific manipulation procedure also shows that 
the trimer-apex tip combined with the single-apex tip has 
potential to achieve single-atom substitutional doping in 
the edge of the cluster and to change its composition, 
which is the motivation of the present work. Usually, the 
edge of the Al nanocluster is modeled by stepped Al (111) 
surface. The extraction and position processes are studied, 
wherein the mechanism is the mechanical interaction force 
acting between the tip apex and surface. An individual 
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atom at the step is extracted first by the tip, and then sin- 
gle Ag or Au dopant is positioned to this site. Based on the 
first-principles simulation, details of the doping process 
are given and its reliability is discussed. 

Methods 

As shown in Figure la, the system we studied is modeled 
by a three-layer Al (111) slab with a step on the top, each 
layer contains 8x7 atoms. The pyramidical Al tip is mim- 
icked by six- or seven-layer atoms mounted on the reverse 
of the slab. In our simulation, two types of tips are consid- 
ered, single-apex tip and trimer-apex one, which are the 
models of sharp and blunt tips in the experiment, respect- 
ively. Besides Al atoms, different tip apex atoms such as 
Ag and Au are taken for doping process. In our survey, 
the tip with different apex atoms can be obtained in the 
experiment [17,18]. As there are six/seven extra layers for 
the tip, our model is convergent with the energy error 
of around 3%. Periodic boundary conditions are imposed 
both parallel and perpendicular to the surface with the 
periodic vectors al, a2, and a3 (see Figure la). By applying 
the periodic boundary condition in Z direction and adjust- 
ing a3, as shown in Figure lb, we can move the tip above 
the stepped surface. The tip is initially placed above the 
step row at a certain height. The distance between the 
tip apex atom and the surface of the lower terrace, 
which is defined as the tip height, is deduced from the Z 
component of the periodic vector a3 (see Figure lb). 

In manipulations, the tip is moved along the X or Z dir- 
ection in a certain step by changing the corresponding 

components of a3 accordingly. Manipulations are simu- 
lated by molecular statics method: after each step, all atoms 
except the ones in the bottom layer of the slab and the top 
layer of the tip are fully relaxed until the forces are smaller 
than 0.01 eV/A. Simulations are based on density func- 
tional theory (DFT) employing the Vienna ab initio simu- 
lation program (VASP) [19]. The exchange-correlation 
potential is described by the generalized gradient approxi- 
mation [20]. Ultrasoft pseudopotentials are used for the 



electron-ion interactions with a cutoff energy of 129 eV 
[21]. The Brillouin zone is sampled with 2 x 4 x 1 k points 
of a Monkhorst-Pack grid. With these parameters, the ob- 
tained lattice parameter of Ag is 4.049 A, which compares 
well with the experimental value of 4.05 A. 

Results 

For the substitutional doping, the first step is extraction of 
surface atom. For this purpose, we consider the trimer- 
apex tip due to its strong attraction to the surface atom 
[11]. Initially, the tip is placed above the manipulated atom 
high enough so that the tip-surface interaction is almost 
negligible, as shown in Figure 2a. Then, we lower down 
the tip step by step. The manipulated atom in the step 
row rises slightly as the tip approaches the surface. When 
the tip height reaches 5.9 A, as shown in Figure 2b, the 
atom is pulled up obviously from the initial site. After that, 
we lift up the tip gradually as shown in Figure 2c to 
Figure 2d; finally, the atom is completely extracted from 
the step site and adsorbed on the tip. During the whole 
process, the tip experiences almost no distortion, which 
indicates that it is stable enough against the atomic inter- 
actions with the surface. In addition, in the extracting 
process, the neighbor atoms of the manipulated atom do 
not show any obvious upward motion, which means that 
the trimer-apex tip can exert effectively attractive force on 
a single atom to make a precise single-atom extraction. 

For understanding the extraction process, as shown in 
Figure 3, we give the total energy varying with the height 
of the manipulated atom relative to the bottom of the 
slab when the tip height is fixed at different heights. 
That is, at a certain tip height, we move the manipulated 
atom down from above in a step of 0.1 A, and at every 
step, the system is relaxed thoroughly. The figure shows 
that at the tip height greater than about 6.3 A, there are 
two local minimum energy wells: one near the surface 
and the other near the tip. When the tip height is lower 
than 6.3 A, the well near the surface disappears grad- 
ually. At 5.9 A, as shown in Figure 3, there is only one 
well near the tip, which means that the manipulated 




Figure 1 Simulation model. The simulation model (a) before and (b) after the periodic boundary condition is applied, in which the tip is 
initially placed above the manipulated atom, (c) The doping process, where the yellow balls represent Al atoms and blue balls represent dopants 
such as Ag or Au. 
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Figure 2 The process of extracting Al atom from the step row 
by the trimer-apex tip. (a) The tip is located upon the 
manipulated atom, (b) Lower down the tip and the manipulated 
atom rises, (c) Lift up the tip gradually, (d) Finally, the atom is 
completely extracted from the step site and adsorbed on the tip. 



atom originally in the step will jump to the well near the 
tip. This is the reason why the manipulated atom obvi- 
ously moves toward the tip when the tip reaches the 
height of 5.9 A. After that, if we lift up the tip, the 
curves in Figure 3 indicate that the manipulated atom 
will stay in the well near the tip. That is, the atom will 



follow the tip and be extracted from the surface, as the 
simulation above shows. From Figure 3, we can also esti- 
mate the reliability of the extraction process; the energy 
curve of 6.1 A shows that the energy barrier for the manip- 
ulated atom escaping from the tip is about 0.25 eV, which 
indicates that the picking up process is robust against the 
disturbances such as thermal diffusion of atoms. 

The next step of substitutional doping is to position a 
dopant atom to the vacancy site where the Al atom is 
extracted. Here, we consider two kinds of dopants: Ag 
and Au atoms. For this purpose, sharp Ag and Au tips 
with single apex atom are considered; such sharp tip can 
be fabricated by electroplating and then annealing, or 
touching a certain metal surface [17,18]. In our simula- 
tion, the sharp Ag tip is modeled by a heterogeneous 
one which contains both Ag and Al atoms, as shown in 
Figure 4. Blue balls indicate the Ag atoms. The apex of 
heterogeneous tip is mimicked by three layers of Ag 
atoms, and our test calculations show that three layers 
of Ag atoms are equivalent to four layers or more. In 
other words, three layers of Ag atoms are sufficient for 
simulation of the sharp Ag tip which is also suitable for 
the Au tip. 

As shown in Figure 4a, the tip is initially placed above 
the vacancy site with the tip height of 8 A at which the 
tip-surface interaction is almost negligible. As the tip ap- 
proaches the surface step by step, the tip apex atom, i.e., 
the dopant atom, relaxes toward the up terrace due to 
the strong attraction. When the tip reaches the height of 
7.1 A, as demonstrated in Figure 4b, the dopant atom 
shows an obvious movement toward the up terrace since 
the attraction is strong enough. At this moment, two 
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Figure 3 Variation of potential energy relative to height of manipulated atom. At different tip heights, the relative potential energy varies 
with the height of the manipulated atom from the Al (111) step surface. 
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Figure 4 The process of positioning Ag dopant to the step site 
by Ag single-apex tip. (a) The tip is located upon the site, (b) As 
the tip approaches the surface, the dopant atom relaxes toward the 
up terrace, (c) Move the tip laterally in the / direction, (d) In the 
end, the dopant atom is released successfully from the tip and 
adsorbed at the step site. 



up-terrace atoms are pulled up slightly and in contact 
with the dopant atom (see Figure 4b). After that, we 
move the tip laterally in the X direction in a step of 
0.2 A at a constant height. As the tip moves forward, as 
shown in Figure 4c, the dopant atom drops gradually be- 
cause of the decreasing vertical attraction from the tip. 
In the end, the dopant atom is released successfully from 
the tip and adsorbed at the step site (see Figure 4d). So far, 
the substitutional doping of the single atom is completed. 

We also give the total energy curves to understand the 
doping process in which the tip laterally moves at a con- 
stant height of 7.1 A. In this case, we fix the tip at different 
lateral distances and move the dopant atom down from 
above in step of 0.1 A, and at every step, the system is re- 
laxed thoroughly. The results, presented in Figure 5, show 
that when the tip stays right upon the vacancy or adsorp- 
tion site, i.e., the lateral distance is 0.0 A, there are two 
local minimum energy wells: one near the surface and the 
other near the tip like the picking up process. The dopant 
atom is still located at the tip because of the energy bar- 
rier. As the tip moves forward along the X direction, the 
right well disappears gradually which means that the at- 
traction from the tip apex is weakened. At the lateral dis- 
tance of 2.4 A, the two wells merge so that the atom 
jumps to the surface. From the curves in Figure 5, it can 
be estimated that the energy barrier for the dopant to es- 
cape from the step site is greater than 0.6 eV, which indi- 
cates that the releasing processes are also reliable even in 
the elevated temperature. 

In order to check the general applicability of our 
substitutional doping method, we next consider the Au 
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Figure 5 Variation of potential energy relative to height of dopant atom. At different lateral distances relative to the vacancy in the X 
direction, the potential energy varies with the height of the dopant atom between the Al (111) surface and the tip. 



Chen et al. Nanoscale Research Letters 2014, 9:235 
http://www.nanoscalereslett.eom/content/9/1/235 



Page 5 of 7 




Figure 6 The process of positioning Au dopant to the Al step 
site by Au single-apex tip. (a) Lower down the tip upon the site, 
(b) Move the tip laterally in the X direction, (c) The Au tip is 
deformed while moving laterally, (d) The dopant atom is 
released successfully. 

V J 



dopant. Similar to the Ag dopant, as shown in Figure 6, 
a single Au atom is also successfully doped into the Al 
stepped surface in the substitutional way. The only dif- 
ference from the case of the Ag dopant is that the Au 
tip is deformed after doping. 



Discussion 

In our doping, both extraction and reposition processes 
only rely on the mechanical interaction force acting be- 
tween the tip apex and the surface. It means that our 
doping scheme, in principle, can be performed with 
STM or AFM. For the STM tip, the electric field is ines- 
sential. Certainly, the specific parameters need to be fur- 
ther confirmed in the experiments. In addition, we find 
that the tip orientation has almost no influence on the 
doping process; as a result, using the tip rotated by 180° 
around the Z axis, we can still achieve the same results. 
The insensitivity to the tip orientation is beneficial to 
the practical experiment. 

We also try other approaches to position the dopant. 
For instance, when the tip reaches 7.1 A, we withdraw 
the tip vertically in the Z direction instead of moving the 
tip laterally in the X direction. For the Ag dopant, it is 
positioned to the vacancy site successfully, as shown in 
Figure 7a,b,c. However, for the Au dopant, it is finally 
placed on the step terrace, as shown in Figure 7d,e,f. 
The difference for Ag and Au can be understood from 
the forces acting on the dopant atom. At the key relax 
step where the dopant atom falls to the surface, we de- 
compose the forces acting on Ag and Au atoms into the 
X and Z directions at every calculation step. The results 
are shown in Figure 8. For the Ag dopant, the compo- 
nent forces have negative peak values and the one in the 
Z direction is greater than that in the X direction, which 
means that the vertical attraction is greater than the lat- 
eral one when the dopant atom is falling. Finally, the Ag 
atom falls into the step site (see Figure 7c). For the Au 
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Figure 8 The forces acting on Ag (a) and Au (b) dopant atom in 
every calculation step. The red curve is the component force in 
the Z direction. The black curve denotes the component force in the 
X direction. 



dopant, however, the component force in the X direction 
has a greater peak value than that in the Z direction. It 
means that the Au dopant tends to drop onto the step 
terrace (see Figure 7f). Though withdrawing the tip ver- 
tically in the Z direction to position the dopant is effect- 
ive for the Ag atom, it lacks general applicability. Also, 
the position details and component forces reveal that it 
is not reliable even in small thermal disturbance (see 
Figures 7 and 8). 

Conclusion 

Based on first-principles simulation, we theoretically inves- 
tigate the substitutional single-atom doping on stepped Al 
(111) surface via atomic manipulation. An effective method 
is proposed in which a trimer-apex tip is adopted to extract 
the surface atom and then a single-apex one is used to pos- 
ition the single dopant atom. In the positioning process, 
the tip moves first in the vertical direction and then in 
a lateral one. Both Ag and Au dopants are successfully 



positioned to the specific site in atomic precision, which 
indicates that the method owns a potential of general 
application. The corresponding energy curves show that 
both extraction and doping processes have a high relia- 
bility against thermal disturbances. Additionally, the ma- 
nipulation processes are insensitive to the tip orientation, 
which is beneficial to the realization of such doping ap- 
proach in practice. 
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